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[0001] This application claims the benefit of the Korean 
Application No. P2003-6730 filed on February 04, 2004, which is 
hereby incorporated by reference. 

BACKGROUND OF THE INVENTION 
Field of the Invention 

[0002] The present invention relates to a display device, and 
more particularly, to an image warping method and apparatus 
thereof, by which image distortion is corrected. 
Discussion of the Related Art 

[0003] Generally, when optical or mechanical distortion is 
caused to such a display device as TV, projector, monitor, and 
the like, image warping uses spatial transformation for 
correcting the distortion. Image warping systems can be 
classified into the following three kinds. 

[0004] 1) Classification according to Transformation Range: 
Global Transformation Method; and Local Transformation Method 

[0005] If coordinates of source image and target image are 
expressed by (u, v) and (x, y) , respectively, the source image is 
represented by the target image, as shown in FIG. 1A, according 
to the global transformation method or the other target image, as 
shown in FIG. IB, according to the local transformation method. 

[0006] Specifically, the global transformation method 

determines spatial transformation positions of all pixels in the 
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image through a polynomial equation expressed by global 
parameters. Hence, the global transformation method has poor 
diversities but enables to provide smooth spatial transformation 
performance without discontinuity using less parameters. 

[0007] On the other hand, the local transformation method is 
performed using a polynomial equation including separate 
parameters for each local area of the image. Hence, post- 
processing is needed since discontinuity may occur at a boundary 
between the local areas. And, the local transformation method 
needs more parameters than the global transformation method since 
separate parameters should be used for each local area. Yet, the 
local transformation method is more advantageous in the 
transformation diversities than the global transformation method. 

[0008] 2) Classification according to Transformation 
Direction: Forward Mapping Method; and Backward Mapping Method 

[0009] A forward mapping method, as shown in FIG. 2, is 
expressed by a transformation relation equation that sets 
coordinates of the source image as independent variables and 
those of the target image as dependent variables, whereas the 
backward mapping method is expressed by another relation equation 
that sets the coordinates of the target image as independent 
variables and those of the source image as dependent variables. 

[0010] In this case, since the forward mapping method maps 

the respective pixels of the source image to the target image, 
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some of the pixels of the target image fail to be mapped (hole 
generation) or are multiply mapped (multiple mapping) . To 
compensate such problems, post -processing such as filtering is 
needed. That's why the backward mapping method is widely used. 

[0011] 3) Classification according to Separability: Separable 
Method; and Non- separable Method 

[0012] Image warping is a sort of two-dimensional spatial 
coordinate transformation, and is a non- separable algorithm in 
horizontal and vertical directions, generally. Yet, many two- 
dimensional transformations can be replaced by continuous linear 
transformation using scanline algorithm (Catmul.l, Edwind and Alvy 
Ray Smith, 3-D Transformations of Images in Scanline Order, 
Computer Graphics, (SIGGRAPH'80 Proceedings), vol. 14, No. 3, 
pp. 279-285, July 1980) . 

[0013] If many two-dimensional transformations can be 
replaced by continuous linear transformation using scanline 
algorithm, it can be regarded as separable in wide sense. 

[0014] FIG. 3 is a block diagram of warping algorithm that is 
horizontally and vertically separable, i.e., scanline algorithm 
proposed by Catmull and Smith. 

[0015] Referring to FIG. 3, a horizontal warping processor 
3 01 receives horizontal scan data and performs warping in a 
horizontal direction to store the result in a memory 302. 



4 



[0016] A vertical warping processor 303 vertically scans to 
read the horizontally warped data stored in the memory 3 02 and 
performs warping in a vertical direction to finally output 
horizontally and vertically warped data. 

[0017] Namely, in case of horizontally/vertically separable 
algorithm, data, as shown in FIG. 3, is processed by horizontal 
and vertical scanning so that a line memory is not needed. 
Moreover, easy data access from memory enables efficient memory 
control . 

[0018] In doing so, processing orders of horizontal/vertical 
warping can be switched. Namely, Catmull and Smith have proposed 
scanline algorithm of forward mapping functions, which is briefly 
explained as follows. 

[0019] First of all, spatial transformation by the forward 
mapping is expressed by Equation 1. 
[0020] [Equation 1] 

[x, y] = T(u, v) = [X(u, v) , Y(u, v) ] , where a function T 
indicates a forward transformation function and functions X 
and Y represent the function T divided by horizontal and 
vertical coordinates , respectively . 
[0021] Hence, once the function T is expressed by T(u, 
v)=F(u)G(v), the function T is separable. In this case, functions 
F and G are called 2 -pass functions. This is because the 
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functions F and G are handled in first and second steps, 
respectively to complete the spatial transformation. 

[0022] However, a general spatial transformation function is 
non- separable . So, the functions F and G become a function of (u, 
v) . Namely, T(u, v) = F(u, v)G(u, v) . 

[0023] For this, Catmull and Smith has proposed the following 
2 -pass algorithm to scanline-process the non-separable function. 

[0024] First of all, if I src , I int/ and I tgt are source image, 
intermediate image, and target image, respectively, 2 -pass 
algorithm can be expressed by the following three steps. 

[002 5] 1 st Step: Assuming that vertical coordinate v of source 
image is constant, a horizontal scanline function can be defined 
as F v (u) = X(u, v) . By performing coordinate transformation 
expressed by Equation 2 in a horizontal direction using the 
mapping function, horizontally warped intermediate image I int is 
made . 

[0026] [Equation 2] 

Isrc(x, V) = I in t(F v (u), V) = I tgt (u, v) 

[0027] Namely, by leaving x v' as assumed constant, data of u 
position of source image is mapped to x position of intermediate 
image . 

[0028] 2 nd Step: The intermediate image prepared by 1 st step 

is represented by (x, v) coordinates. In doing so, since the 

intermediate image expressed by (u, v) coordinates are needed for 
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horizontal processing, an auxiliary function H x (v) is driven by 
adjusting x = X(u, v) of Equation 1 for *u' by leaving *x' as 
constant. Namely, as u = H x (v), it is represented by a function 
of x v' . In this case, the auxiliary function H x (v) is usually not 
expressed as a closed form. In such a case, such a numerical 
method as Newton -Raphson iteration method is needed to solve the 
equation. 

[0029] 3 rd Step: A vertical scanline function is defined as 
follows using the auxiliary function. As G x (v) = Y(H x (v), v) , a 
function of x v' only is prepared. So, warping can be executed in 
a vertical direction. Namely, coordinate transformation expressed 
by Equation 3 is executed by scanning the intermediate image in a 
vertical direction using the mapping function of the variable v, 
i.e., G x (v) = Y(H x (v), v) , thereby providing the 
horizontally/vertically warped target image I tgt . 

[0030] [Equation 3] 

Itgt(x, y) = I tgt (x, G X (V)) = Iint(x, v) 

[0031] In this case, the most difficult work in implementing 
the scanline algorithm is to seek the auxiliary function of the 
2 nd step since it is generally unable to find an auxiliary 
function expression of closed form. 

[0032] Hence, in U.S. Patent No. 5,204,944 (George Wolberg, 

Terrance E. Boult, Separable Image Warping Methods and Systems 

Using Spatial Lookup Tables) , disclosed is a method of 
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implementing the above -explained scanline algorithm for the local 
transformation method and the forward transformation method. In 
this case, input coordinates are simultaneously re-sampled 
together with image data using a lookup table to solve the 
problem of finding the auxiliary function. 

[0033] However, the above method needs excessive hardware for 
re-sampling coordinates. Moreover, as mentioned in the foregoing 
description of the forward transformation method, mapping failure 
(hole generation) or multiple mapping of the' target image pixels 
may occur. Hence, the forward transformation method needs 
separate post -processing to raise algorithm complexity, thereby 
increasing implementation costs. 

SUMMARY OF THE INVENTION 

[0034] Accordingly, the present invention is directed to an 
image warping method and apparatus thereof that substantially 
obviates one or more problems due to limitations and 
disadvantages of the related art. 

[003 5] An object of the present invention is to provide an 
image warping method and apparatus thereof, by which simplified 
scanline algorithm is implemented by a backward transformation 
method with minimized implementation costs and which enables to 
correct image distortion of a display device such as projection 
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TV, projector, monitor, and the like due to optical or mechanical 
distortion. 

[0036] Additional advantages, objects, and features of the 
invention will be set forth in part in the description which 
follows and in part will become apparent to those having ordinary 
skill in the art upon examination of the following or may be 
learned from practice of the invention. The objectives and other 
advantages of the invention may be realized and attained by the 
structure particularly pointed out in the written description and 
claims hereof as well as the appended drawings. 

[0037] To achieve these objects and other advantages and in 
accordance with the purpose of the invention, as embodied and 
broadly described herein, an image warping method according to 
the present invention includes a step (a) of, if coordinates of 
source and target images are defined as (u, v) and (x, y) , 
respectively, driving an auxiliary function by finding a solution 
of the coordinate y of the target image by leaving the coordinate 
v of the source image as constant, a step (b) of preparing a 
horizontally warped intermediate image by applying the auxiliary 
function to a first backward mapping function u = U(x, y) , and a 
step (c) of preparing a horizontally/vertically warped target 
image by applying the horizontally warped intermediate image to a 
second backward mapping function v = V(x, y) . 
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[0038] In this case, the step (b) includes a step (d) of 
finding the coordinate u of the source image by receiving to 
apply a value of the coordinate x of the target image, polynomial 
coef f icient (s) of the first backward mapping function, and the 
auxiliary function to the first backward mapping function and a 
step (e) of preparing the horizontally warped intermediate image 
by interpolating data of the coordinate u found in the step (d) . 

[0039] And, the step (c) includes a step (f) of applying the 
second backward mapping function to the intermediate image, a 
step (g) of finding the coordinate v of the source image by 
receiving to apply values of the coordinates x and y of the 
target image, polynomial coef f icient (s) of the first backward 
mapping function, and a result applied in the step (f) to the 
second backward mapping function, and a step (h) of preparing a 
horizontally/vertically warped target image by interpolating data 
of the coordinate v found in the step (g) . 

In another aspect of the present invention, an image warping 

method includes a step (a) of, if coordinates of source and 

target images are defined as (u, v) and (x, ' y) , respectively, 

driving an auxiliary function (y = H v (x) ) from a backward 

mapping function v = V(x, y) by finding a solution of the 

coordinate y of the target image by leaving the coordinate v of 

the source image as constant, a step (b) of preparing a 

horizontally warped intermediate image by applying the auxiliary 
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function (y = H v (x) ) to a backward mapping function u = U(x, y) , 
and a step (c) of preparing a horizontally/vertically warped 
target image by applying the horizontally warped intermediate 
image to the backward mapping function v = V(x, y) . 

[0040] In another aspect of the present invention, an image 
mapping apparatus includes a horizontal warping processing unit 
providing a horizontally warped intermediate image, if 
coordinates of source and target images are defined as (u, v) and 
(x, y) , respectively, by receiving a value of the coordinate x of 
the horizontally scanned target image and coefficients b 0 o~b 2 i of 
a polynomial, by finding a solution of the coordinate y of the 
target image by leaving v as constant to drive an auxiliary 
function (y = H v (x)), and by applying the auxiliary function (y = 
H v (x)) to a backward mapping function u = U(x, y) , a memory 
storing the horizontally warped intermediate image of the 
horizontal warping processing unit, and a vertical warping 
processing unit providing a horizontally/vertically warped target 
image by scanning the horizontally warped intermediate image 
stored in the memory in a vertical direction and by applying the 
scanned image to a backward mapping function v = V(x, y) . 

[0041] In this case, the horizontal warping processing unit 

includes a first auxiliary function computing unit driving the 

auxiliary function (i.e., Ay 2 +By+C=0) by receiving the value of 

the coordinate x of the horizontally scanned target image and the 
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coefficients b 0 o~b 2 i of the polynomial and by adjusting backward 

mapping function for y by leaving v as constant, a second 

auxiliary function computing unit finding a solution (y = H v (x)) 

for the auxiliary function, a u-coordinate computing unit finding 

the coordinate u of the source image by receiving the coordinate 

x of the target image, coefficients a 0 o~a 2 i of the polynomial, and 

a value of y for the auxiliary function, an address and 

interpolation coefficient detecting unit outputting an integer 

part u int of the coordinate u as an address assigning a data-read 

position in the memory and a fraction part ( a = u - u int ) as an 

interpolation coefficient, and an interpolation unit 

interpolating data I S rc(Ui nt/ v) of the source image outputted from 

the memory with the interpolation coefficient a to output the 

interpolated data as the intermediate image Ii n t(x, v) . 

[0042] And, the vertical warping processing unit includes a 

v- coordinate computing unit finding the coordinate v of the 

source image by scanning the intermediate image stored in the 

memory and by receiving x and y of the target image and the 

coefficients b 00 ~b 2 i of the polynomial, .an address and 

interpolation coefficient detecting unit outputting an integer 

part v int of the coordinate v as an address assigning a data-read 

position in the memory and a fraction part a {a = v - v int ) as an 

interpolation coefficient, and an interpolation unit outputting 

the target image I tg t(x, y) by interpolating data Ii n t(x, v int ) of 

12 



the intermediate image outputted from the memory with the 
interpolation coefficient a outputted from the address and 
interpolation coefficient detecting unit. 

[0043] It is to be understood that both the foregoing general 
description and the following detailed description of the present 
invention are exemplary and explanatory and are intended to 
provide further explanation of the invention as claimed. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0044] The accompanying drawings, which are included to 
provide a further understanding of the invention and are 
incorporated in and constitute a part of this application, 
illustrate embodiment (s) of the invention and together with the 
description serve to explain the principle of the invention. In 
the drawings : 

[0045] FIG. 1A and FIG. IB are diagrams of global and local 
transformation methods of image warping, respectively; 

[0046] FIG. 2 is a diagram of forward and backward mapping 
methods of image warping; 

[0047] FIG. 3 is a block diagram of warping algorithm that is 
horizontally and vertically separable; 

[0048] FIGs. 4A to 4H are diagrams of distortion types 
existing on a general display device; 



13 



[0049] FIG. 5 is a block diagram of a horizontal warping 
processor according to the present invention; 

[0050] FIG. 6 is a block diagram of a vertical warping 
processor according to the present invention; and 

[0051] FIG. 7 is a diagram of a bilinear interpolation method 
according to the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 
[0052] Reference will now be made in detail to the preferred 
embodiments of the present invention, examples of which are 
illustrated in the accompanying drawings. Wherever possible, the 
same reference numbers will be used throughout the drawings to 
refer to the same or like parts. 

[0053] Geometrical spatial transformation is generally needed 
to correct image distortion caused to an image display device by 
optical or mechanical factors. In this case, when coordinates of 
source and target images are expressed by (u, v) and (x, y) , 
respectively, a backward mapping function used for spatial 
transformation has such a polynomial form as Equation 4. 
[0054] [Equation 4] 
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u = U(x,y) = f d ^a ij x i y j 

i=0 7=0 

j. , . , , , where a i and b { are coefficients of 

/V N-i 1 1 

1=0 y=0 

polynomial, respectively and N indicates an order of 
polynomial . 

[0055] In this case, there exist more distortion types as the 
order of the polynomial increases. Yet, as the coefficient of the 
polynomial increases, algorithm complexity, and implementation 
costs are raised. 

[0056] And, the distortion types appearing on the display 
device are shown in FIGs . 4A to 4H. In this case, correctable 
distortions are explained according to the order of polynomial as 
follows . 

[0057] When the polynomial order is x l', there are shifting 
(FIG. 4A) , scaling (FIG. 4B) , horizontal skew (FIG. 4C) , vertical 
skew (FIG. 4D) , and tilt (FIG. 4B) . 

[0058] When the polynomial order is x 2', there is keystone 
(FIG. 4E) . 

[0059] When the polynomial order is '3', there are pincushion 
(FIG. 4G) and barrel (FIG. 4H) . 

[0060] Hence, in order to correct the distortion types shown 
in FIGs. 4A to 4H, at least cubic polynomial should be used. 
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[0061] Meanwhile, scanline algorithm of a backward mapping 
function proposed by the present invention is executed by the 
following three steps. 

[0062] 1 st Step: An auxiliary function H v (x) is driven from 
v=V(x, y) of Equation 4 by finding a solution for a vertical 
coordinate % y' of a target image by leaving a vertical coordinate 
*x' of a source image as constant. Namely, y=H v (x). 

[0063] 2 nd Step: If the auxiliary function H v (x) found by 1 st 
Step is applied to the first function of Equation 4, u=U(x, 
H v (x)), which is represented by a function of x x' only. Using 
this mapping function, a horizontally warped intermediate image 
lint is provided by Equation 5. 

[0064] [Equation 5] 

IintU, y) = I src (U(x, H V (X)), V) = I sr c(u, v) 

[0065] 3 rd Step: A horizontally/vertically warped image is 
provided by Equation 6 using the second function, v=V(x, y) , of 
Equation 4 applied to the horizontally warped image of Equation 5. 

[0066] [Equation 6] 

Itgt(x, y) = Iint(x, V(X, y) ) = I in t(x, v) 

[0067] Namely, after a position x u' of the source image has 

been found using the value of x x' of the target image, data of 

the position x u' of the source image is mapped to a position x x' 

of the target image. After a position x v' of the source image has 
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been found using the values of *x' and x y' of the target image, 
data of the position W is brought to be mapped to a position 
% y' of the target image. Hence, horizontally/vertically warped 
data can be attained. In doing so, the sequence of the horizontal 
and vertical warping can be switched. 

[0068] An image warping method, which is implemented from the 
scanline algorithm proposed by the present invention and the 
global mapping function expressed by Equation 1, is explained in 
detail as follows. 

[0069] First Embodiment 

[0070] The cubic polynomial developed from Equation 4 is 
represented by Equation 7. 
[0071] [Equation 7] 

u = U(x,y) = a 00 + a ox y + a 02 y 2 + a 03 y 2 + a XQ x + a x x xy + a x2 xy 2 + a 20 x 2 + a 2X x 2 y + a 30 x* 
v = V(x,y) = b 00 + b Qx y + b 02 y 2 +b 03 y 3 + b l0 x + b u xy + b X2 xy 2 +b 2Q x 2 + b 2x x 2 y + b 3Q x> 

[0072] If cubic terms in Equation 7, which are unnecessary 
for correcting the distortion types shown in. FIGs . 4A to 4H, are 
removed, a mapping function shortened as Equation 8 is attained. 

[0073] [Equation 8] 

u = U(x, y) = a 00 + a ox y + a 02 y 2 + a xo x + a u xy + a n xy 2 + a 20 x 2 + a 2X x 2 y 
v = V(x, y) = bn+ b ox y + b Q2 y 2 + b xo x + b xx xy + b n xy 2 + b 20 x 2 + & 2l x 2 j/ 

[0074] In order to calculate an auxiliary function, if the 

second one of Equation 8 is adjusted for y by leaving x v' as 

constant, a quadratic function is represented by Equation 9. 
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[0075] [Equation 9] 

Ay 2 + By + C = 0 , where A=b 02 + b n x , B=b 0l +b u x + b 2l x 2 , and 

C=b 00 +b w x + b 20 x 2 - v . 

[0076] Hence, from root formula, % y' of Equation 9 can be 
driven as Equation 10. 
[0077] [Equation 10] 

y-H v (x)- 

[0078] In this case, there may exist three kinds of roots in 
Equation 10. And, a processing method should vary according to 
each case. 

[0079] First of all, if there exist two real roots ( B 2 > 4 AC) , 



+ -B + Jb 2 -4AC -B-Jb 2 -4AC 

one of two rear roots, v = - and v = , 

' 2A y 2A 

is arbitrarily selected to use. 

[0080] And, in case that there exists one real root (B 2 = 4 AC) , 
— B 

the root is y = . 

2A 

[0081] Moreover, if there exist a pair of imaginary roots 



t D 2 * ^ + -B + iJB 2 -4AC 
{ B < 4 AC) , a pair of the imaginary roots are y = 



2A 



_ -B-iJB 2 -4AC 

and y = . in this case, if B < 4 AC, since 

2A 

coordinates having imaginary values substantially fail to exist, 
the imaginary terms of the equation are ignored to provide the 
same root of the second case that there exists one real root. 
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[0082] After the auxiliary function, y = H v (x), has been 
calculated, spatial transformations are horizontally and 
vertically executed using Equations 5 and 6, respectively. In 
doing so, since the coordinates mapped by the transformation 
function are not located at the pixel sample of the source image 
in general, a pixel value of the mapped coordinates is calculated 
by interpolation using neighbor pixels. 

[0083] Specifically, assuming that a center of image is set 
as an origin of coordinates and that sizes of source and target 

images are set to W src x H src and W tgt x H tgt , respectively, input 

coordinates v and x of a horizontal warping processor are 

—H H —W W 

integers having ranges of [ — , — — ] and [ — , — — ], 

2 2 "22 

respectively. And, a coordinate x u' calculated by the horizontal 
warping processor becomes a real number. In doing so, an integer 
part u m is used as an address for assigning a location of data 
to read from a memory and a fraction part a is used as an 
interpolation coefficient. 

[0084] A first auxiliary function computing unit 501, as 
shown in FIG. 5, of the horizontal warping processor receives x x' 
of the horizontally scanned target image and coefficients boo~b 2 i, 
and adjusts the polynomial for s y' like Equation 9 by leaving *v' 
as constant to express a quadratic function (i.e., 
Ay 2 + By + C = 0). And, a second auxiliary function computing unit 
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502 finds a solution of x y' , i.e., auxiliary function, like 
Equation 10 (y = H v (x)) using the root formula to output to a u- 
coordinate computing unit 503. 

[0085] The u-coordinate computing unit 503 applies the 
auxiliary function to the first function of Equation 4 to find a 
coordinate x u' of the source image by receiving x x' , a 0 o~a 2 i, and 
x y' found by Equation 10, and then outputs the coordinate x u' to 
an address and interpolation coefficient detecting unit 504. 

[0086] The address and interpolation coefficient detecting 
unit 504 outputs the integer part w int of the coordinate x u' as an 
address for assigning a location of data to read to a memory 505 
and the fraction part {a = u - u int ) as an interpolation 
coefficient to an interpolation unit 506. • - 

[0087] The memory 505 outputs data I S rc(u in t/ v) of the source 
image stored in the inputted address addr to the interpolation 
unit 506. And, the interpolation unit 506 interpolates the data 
Isrc(u in t, v) of the source image outputted from the memory 505 
with the interpolation coefficient a outputted from the address 
and interpolation coefficient detecting unit 504, thereby 
outputting the intermediate image I int (x, v) like Equation 5. 

[0088] In this case, since the coordinates mapped by the 
transformation function of Equation 5 is not located at the pixel 
sample u of the source image in general, the interpolation unit 
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506 calculates the pixel value of the mapped coordinates by 
interpolation using neighbor pixels. 

[0089] FIG. 7 is a diagram of a method using bilinear 
interpolation. 

[0090] Namely, bilinear interpolation in FIG. 7 can be 
represented by Equation 11. 
[0091] [Equation 11] 

/ int (x, v) = (1 - a) I src (w int , v) + a I src (w int - 1, v) 

[0092] If the horizontal warping image processor in FIG. 5 is 
firstly operated, the horizontally warped intermediate image is 
stored in the memory. Thereafter, the intermediate image stored 
in the memory is scanned in a vertical direction and is then 
applied to the backward mapping function to provide the 
horizontally/vertically warped target image, finally. 

[0093] Meanwhile, referring to a vertical warping image 
processor of FIG. 6, a v-coordinate computing unit 601 scans the 
intermediate image stored in the memory in a vertical direction, 
finds the coordinate x v' of the source image by receiving to 
apply x and y of the target image and the coefficients b 0 o-b 2 i of 
the polynomial to Equation 7, and then outputs the found 
coordinate *v' to an address and interpolation coefficient 
detecting unit 602. 

[0094] 
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[0095] The address and interpolation coefficient detecting 
unit 602 outputs the integer part v int of the coordinate x v' as an 
address for assigning a location of data to read to a memory 603 
and the fraction part a ( a = v - v int ) as an interpolation 
coefficient to an interpolation unit 604. 

[0096] The memory 603 outputs data Iint(x, v int ) of the 
intermediate image stored in the inputted address addr to the 
interpolation unit 604. And, the interpolation unit 604 
interpolates the data Ii n t(x, v int ) of the intermediate image 
outputted from the memory 603 with the interpolation coefficient 
a outputted from the address and interpolation coefficient 
detecting unit 603, thereby outputting the target image I tg t(x, y) 
like Equation 6. 

[0097] Likewise, since the coordinates mapped by the 
transformation function of Equation 6 is not located at the pixel 
sample v of the source image in general, the interpolation unit 
604 calculates the pixel value of the mapped coordinates by 
interpolation using neighbor pixels. 

[0098] Second Embodiment 

[0099] In the first embodiment of the present invention, the 

part for computing the auxiliary function needs relatively 

excessive calculation load and hardware. By adopting small 

approximation, such calculation load and hardware can be reduced 

without degrading warping performance. 
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[00100] Namely, in most cases of the quadratic function of 
Equation 9, X A' is much smaller, than X B' or X C . Using such a 
fact, Equation 9 can be approximated to a linear function of 
Equation 12. 

[00101] [Equation 12] 

By + C = 0 , where B=b 0{ +b u x + b 2l x 2 and C=b 00 + b w x + b 20 x 2 -v. 

[00102] From the root formula, x y' of Equation 12 can be 
simply found as Equation 13. 
[00103] [Equation 13] 

y = H v (x) = ^ 

[00104] After the auxiliary function y = H v (x) has been 
calculated, horizontal and vertical transformations are executed 
using Equations 5 and 6. 

[00105] Accordingly, an image warping method and apparatus 
thereof according to the present invention enables to implement 
the simplified scanline algorithm by the backward transformation 
method with minimized implementation costs and to correct the 
image distortion, which is caused by optical or mechanical 
distortion, of a display device such as projection TV, projector, 
monitor, and the like. 

[00106] Namely, the present invention adopts the backward 
mapping method to avoid pixels of non-mapping or multiple mapping, 
and uses the global transformation method to enable the smooth 
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spatial transformation without discontinuity for the entire image 
with less parameters. Therefore, the present invention needs no 
additional post-processing. 

[00107] And, by adopting the scanline algorithm, the present 
invention enables the efficient memory access as well as 
simplified circuit configuration and cost reduction in aspect of 
hardware implementation. Therefore, the present invention is very 
competitive in costs and performance when applied to display 
devices such as projection TV, projector, monitor, etc. to which 
the correction of image distortion is essential. 

[00108] It will be apparent to those skilled in the art that 
various modifications and variations can be made in the present 
invention. Thus, it is intended that the present invention covers 
the modifications and variations of this invention provided they 
come within the scope of the appended claims and their 
equivalents . 
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